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A study is made of some properties of this interface in the SU(3) pure gauge theory in 2+1 dimensions. At high 
temperatures, the interface tension is measured and shows agreement with the perturbative prediction. Near the 
critical temperature, the behaviour of the interface is examined, and its fluctuations compared to a scalar field 
theory model. 
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1. Preliminaries 

The pure glue sector of QCD, SU(3) gauge 
theory, possesses several phases at finite tem- 
peratures (T). At low T, it exists in a single 
disordered phase where the vacuum is symmet- 
ric under the centre of the group, Z{2>), and 
the colour charge is confined; at high T, there 
are three distinct ordered, colour-non-confining, 
Z(3)-breaking phases, each degenerate vacuum 
corresponding to an element of Z{2>). At some 
critical temperature, T c , there is a phase transi- 
tion between the two regimes. Thus, two types of 
phase interface can occur in the theory: one be- 
tween ordered and disordered phases, stable only 
at T c , and one between two ordered phases, sta- 
ble for T > T c . The latter is often referred to as 
an "order-order" or "Z(3)" interface. 

In the Euclidean formalism, the "time" dimen- 
sion runs from to l/T, and an order parameter 
for the phases is the Polyakov line, a time-ordered 
Wilson loop which wraps around the time bound- 
ary for a fixed spatial location f: 
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The expectation value of this Polyakov line gives 
the self-energy of a single quark in the gluon 
medium: < L >= e~ F l T . Confinement implies 
F = oo so that < L >= in the disordered phase. 
For each ordered phase, < L > is equal to a differ- 
ent element of Z(3)[l]. Unlike topologically triv- 
ial Wilson loops, then, the Polyakov line is aware 
of the phase structure of an SU (3) system and 



can map out the structure of any interface. 

It is widely believed, and has been shown by 
simulation [2], that the interfaces of SU{2>) show 
"critical wetting" , meaning that a Z{2>) interface 
should be thought of as two order-disorder in- 
terfaces stuck back-to-back, leaving a thin slice 
of disorder within. Possible cosmological impli- 
cations have been discussed in the past [3], with 
Z{2>) interfaces from the hot quark-gluon plasma 
splitting apart, as the early universe cooled to 
T c , to spread confined matter outwards. How- 
ever, questions have been raised as to whether 
these essentially Euclidean objects can have any 
significance in the real universe, especially as the 
high-T order leads to unphysical thermodynamic 
properties for the interfaces[4] in the presence of 
fermions, and it has been argued that there may 
be only one true Z{2>) phase even at high T[5]. 

2. The Interface Tension at High T 

The calculation of the interface tension, a, is 
an instanton problem, in which the Polyakov line 
interpolates between different Z{2>) vacua. An ef- 
fective one-dimensional theory describes the in- 
terface profile, with the effective action for the 
Polyakov line consisting of a classical kinetic 
part and a potential term formed by integrat- 
ing out one-loop quantum fluctuations about a 
constant background field. The potential incor- 
porates the Z{2>) symmetry, having minima when 
L G Z(3)[l]. The instanton calculation has been 
performed[6], predicting a behaviour in 2+1 Eu- 
clidean dimensions of the form 
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where «o is a constant predicted to be 8.33 for 
SU(3) in the continuum, in the limit as T —> 
oo. On the lattice, this prediction is modified 
to 9.82. It has also been argued[5] that higher- 
order infrared divergences invalidate this calcula- 
tion, leaving only one true Z{2>) phase and, con- 
sequently, an interface tension of zero. 

In order to test this prediction, we simulate 
SU{2>) pure gauge theory in 2+1 dimensions on a 
cubic Euclidean lattice of dimension L z xL x xL t , 
L z > L x Lf. The familiar Wilson plaquette ac- 
tion is used, with the variables defined on the lat- 
tice links being SU(3) matrices. A "twist" [2] en- 
sures the appearance of a physical interface: a set 
of plaquettes pierced by a skewer in the x direc- 
tion has its contribution to the action multiplied 
by an element of Z{2>), representing an effective 
change of variables; periodic boundary conditions 
then require the field to interpolate between the 
two resulting vacua elsewhere in the z direction. 

Simulating with and without the twist and 
measuring the average action in each case gives 
us the difference AS, from which we obtain AF, 
the free energy of the interface. With lattice spac- 
ing a, the interface "area" is A = L x a and the 
temperature T = l/(_L t a), so that 

1 a c 8 ( aA \ L * (*\ 

where the second equality uses (2) and introduces 
the lattice parameter (3 = 6/(ag 2 ). 

Our simulations were performed for various 
16.0 < j3 < 200.0 with L t = 2. The spatial di- 
mensions were scaled with the interface width, 
~ \/]3; for (3 = 200, we used L z = 96 and 
L x = 32. To estimate ao, it is necessary to ex- 
trapolate to [3 = oo. Unfortunately, the exact 
form of the perturbative correction to a in (2) is 
not known. In 3+1 dimensions, a 1/ [3 dependence 
is seen, but there is good reason to expect addi- 
tional factors of In [3 from the infrared divergences 
in 2+1 dimensions, as seen in the Debye mass[7]. 
In fig. 1, we perform the extrapolation using a 
best fit to a power of 1//3, which we measure to 
be 0.85(10). The extrapolated value obtained is 
a = 9.91(1)(14), the second error being asso- 
ciated with the power of 1//3. This value is in 
good agreement with the instanton calculation, 



suggesting that the Z{2>) phases are indeed dis- 
tinct at high T, and must be taken into account 
in calculations from the Euclidean path integral. 
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Figure 1. A linear extrapolation of «o against 
1/13° 85 to give the j3 = oo limit. 

3. Qualitative Interface Behaviour Near T c 

At high T ^> T c , the two Z{2>) phases intro- 
duced on the lattice by use of the twist are well 
defined and separated by a rigid interface; fluctu- 
ations in phase are very small. However, at and 
below T c , these Z{2>) phases cease to exist, along 
with any Z{2>) interface, being replaced by the 
single disordered phase with symmetry restored. 
This raises the question: how does the interface 
behave between these two regimes? 

As T decreases towards T c , a corresponding de- 
crease is seen in the energy penalty suppressing 
phase fluctuations. One also expects the expecta- 
tion values of Polyakov lines to converge on zero 
at the (second-order) phase transition, this being 
the value seen in the disordered phase. The com- 
bination of these factors means that the "height" 
of the Z{2>) interface (the difference in the real 
part of the Polyakov line as one moves across it) 
decreases just as more and larger bubbles of dif- 
ferent phase appear within the two main Z{2>) do- 
mains on the lattice. The interface thus becomes 
increasingly difficult to discern as one approaches 
T c . For our system, it is hard to see the interface 
on a plot of Polyakov line values below [3 ps 9.00, 
against a critical value of f3 c ps 8. 175 [8]. This, 
of course, is precisely the temperature range of 



interest. 

The first objective must be to uncover the qual- 
itative behaviour of the interface as T drops. 
Several possibilities present themselves, given its 
properties at high T and at T c : the interface could 
maintain its high-T rigid structure with mini- 
mal "area" , its "height" shrinking steadily until 
it vanishes at T c ; it could spread out longitudi- 
nally as correlation lengths on the lattice diverge 
with decreasing T, revealing more and more dis- 
ordered phase and becoming less and less mean- 
ingful as a physical object; it could be washed out 
by the formation of many more interfaces along 
with the large bubbles of phase, colliding with 
each other and leaving disordered phase as the 
energy penalty drops; or it could maintain a rel- 
atively narrow width, but suffer increasingly vio- 
lent transverse fluctuations which diverge towards 
T c , where it would break apart to leave only dis- 
order. 

To follow the interface into the regime near T c , 
one needs a way to identify it amidst the sea of 
phase fluctuations. One such way is to screen out 
the highest-frequency fluctuations, for instance 
by a straightforward box-car average over neigh- 
bouring Polyakov lines. This produces a spatial 
phase profile where the interface can be seen much 
more clearly, at the cost of some loss of informa- 
tion about its shape. The interface becomes vis- 
ible as a localised, narrow object right down to 
T e . 

An improvement on this method is to contour 
the Polyakov profiles, picking various "heights" 
between the two extremes at opposite longitudi- 
nal ends of the lattice. Assuming no other inter- 
faces are produced on the lattice, only one con- 
tour at each level will wrap across the lattice, hav- 
ing unit "winding number" ; mere phase bubbles 
will be represented by "topologically trivial" con- 
tours. This produces a detailed picture of the 
phase structure of the interface at any stage of 
the simulation. Our results reveal that the inter- 
face does indeed remain thin and localised, but 
with increasingly large fluctuations in shape as 
T drops, as illustrated in the snapshot of fig. 
2. This matches the fourth of the possible be- 
haviours listed previously. 
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Figure 2. A contour snapshot of Polyakov lines 
(real parts) using three contour levels for L x = 
36, = 72 and /? = 8.5. 

4. Quantitative Analysis 

Having ascertained the qualitative behaviour of 
the interface to be that of a fluctuating string 
whose fluctuations increase as T —> T c , the next 
task is to gain a better quantitative understand- 
ing. One might suspect that the fluctuations 
would diverge at T c itself, dissolving the interface 
and leaving only disorder, as desired. To test this, 
we follow various fluctuation moments: if <j){x) is 
the displacement from average of an interface con- 
tour at transverse location x, the n'th moment is 
defined to be <(f)(x) n >. We track the first six mo- 
ments, taking a separate average for each trans- 
verse position. Since the interface is translation 
invariant, with no preferred longitudinal position, 
one expects the moments to be roughly constant 
across the lattice after a large number of sweeps, 
and this is indeed what we see, with the even 
moments also found to be much larger than the 
odd. Thus, we can further average the moments 
across the lattice, and by producing such aver- 
age moments for various temperatures near T c , 
we see that the even moments do indeed diverge 
as T —> T c . 

The behaviour of the interface, rigid at high T 
but with increasing fluctuations as T drops, sug- 
gests a description in terms of a scalar field theory 
with a simple interaction to represent production 
and destruction of bubbles of phase at the in- 
terface. To maintain translation invariance, and 
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Figure 3. A log plot of the ;c-averaged, Wick- 
subtracted, n'th (even) moments as a function of 
L x , for /3 = 8.5. 

produce the domination of the even moments, we 
suggest a Lagrangian of the following form: 

C = i 7 (^) 2 " U^) 4 - (4) 

Here, j and A are unknown functions of (/? — f3 c ). 
We then equate the n'th fluctuation moment, af- 
ter Wick subtraction, to the n-point connected 
vacuum correlation function of this Lagrangian. 
As desired, the odd moments are predicted to be 
zero, since the interaction is of even order. A test 
of the applicability of the Lagrangian is provided 
by the dependence of the moments on the lattice 
width (the finite size of the string). Correlation 
functions predict dependences of L x , L x and L x 
for the 2nd, 4th and 6th moments respectively, 
and this is roughly what we see in results such as 
fig. 3, with slopes of 1.16(1), 2.02(5) and 1.8(4). 

Having satisfied ourselves that the scalar field 
theory model provides a reasonable description of 
the interface near its breakdown, we can measure 
the moments at different T to find the functions 
7 and A. For 8.25 < [3 < 9.00, we obtain results 
such as fig. 4. From the gradients, fitted to be 
-0.53(1), -1.59(3) and -2.9(4), we estimate that 
7 ~ (J3 - /3 C )° 53 (!) and A ~ (j3 - /? c )° 53 ( 4 ). 

5. Conclusions 

We have measured the tension of a Z{2>) in- 
terface in the T = oo limit, and found excellent 
agreement with theory, bolstering the argument 
for the existence of Z(3)-breaking phases in the 
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Figure 4. A log plot of the ;c-averaged, Wick- 
subtracted, n'th (even) moments as a function of 
(/?-&), for L x = 36. 

Euclidean path integral. We have also studied 
the behaviour of this interface in a very differ- 
ent regime, near to its collapse. Measuring very 
close to the critical temperature, we have found 
its qualitative behaviour to be that of a string 
with divergent fluctuations at the critical point, 
these being well described by a simple interacting 
scalar field theory. Further analysis in continuing, 
and a similar study in 3+1 dimensions may well 
be rewarding. 
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